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Abstract 

Determining contributions by sub-portfolios or single exposures to portfolio- wide economic capital 
for credit risk is an important risk measurement task. Often economic capital is measured as Value- 
at-Risk (VaR) of the portfolio loss distribution. For many of the credit portfolio risk models used in 
practice, the VaR contributions then have to be estimated from Monte Carlo samples. In the context 
of a partly continuous loss distribution (i.e. continuous except for a positive point mass on zero), we 
investigate how to combine kernel estimation methods with importance sampling to achieve more 
efficient (i.e. less volatile) estimation of VaR contributions. 

1 Introduction 

In many financial institutions, there is a well esta blished practice of measuring the risk of their portfolios 



in terms of economic capital (cf., e.g. iDevl . 120041 ). Measuring portfolio-wide economic capital, however, 
is only the first step towards active, portfolio-oriented risk management. For purposes like identification 
of risk concentrations, risk-sensitive pricing or portfolio optimisation it is also necessary to decompose 
port folio- wide econom ic capital into a sum of risk contributions by sub-portfolios or single exposures (see, 
e.g- lLitte rmanl . 119961 ). 

While already calculating or estimating economi c capital is non-trivial in g eneral, determining risk de- 



compositions is even more demanding (see, e.g.. lYamai and Yoshibal . 120021 ) . For most of the economic 



capital models used in practice, no closed-form solutions for risk contributions are available^- As a conse- 
quence, for such models there is a need for simulation to create samples from which to estimate economic 
capital as well as risk contributions. These estimations often involve evaluations of very far tails of the 
risk return distributions, causing high variability of the estimates. Various var iance reduction techniques 



have been proposed , one of the more popular bein g importance sampling (see iGlasserman and Lit I200E 
Merino and Nvfelerl . |2004| ; iKalkbrener et al. , 2004, for its application to credit risk). IGlasserman! (|2005 



suggests a two-step importance sampling approach to the estimation of contributions to Value-at-Risk 
(VaR), the most popular metric underlying economic capital methodologies. 

IGlasserman 's approach, however, does apply to discrete loss distributions only (i.e. to distributions such 



that each potential loss value has a positive probability to be assumed), which means a significant 
restriction as loss distributions based on continuous loss given default rate distributions seem to be more 
realistic in practical applications. According to Theorem 13.31 below such loss distributions have the 
property that each single potential positive loss has probability zero to be assumed but there is a positive 
probability of not observing any loss. The probability of not observing losses is significantly positive in 
particular for small and medium portfolio sizes (with 200 names or less) which occur, for instance, in 
typical securitisation deals. 
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To solve the problem of determining VaR contributions whe n positive losses have a continuous distribu- 
tiorfl in thi s paper we follow, wh e re pos sible, the path used bv lGourieroux et al.l (|2000 ), lEpperlein and Smillid 
( 20061 ). and Gourieroux and Liu ( 20061 ) who apply kernel estimation methods for estimating VaR contri- 
butions and contributions to spectral risk measures in a ma rket risk context with co ntinuous distributions. 
Kernel estimators are a well-established concept (see, e.g., Pagan and Ullab . 19991 ) to deal with the issue 
of estimating non-elementary expectations as they occur in the context of capital allocation. Due to the 
rare event issue characteristic for credit risk, entailing rather volatile estimates when standard Monte 
Carlo simulation is used, we combine the kernel estimation technique with importance sampling (shifting 
the means of the systematic factors to be more specific) for reducing estimation variance. The paper is 
organized as follows: 



Scction[2]contains a review, in the necessary details, of the capital allocation problem and the specific 
issues with the estimation of risk contributions when the loss distribution is partly continuous. 

Section [3J provides a brief review of kernel estimators for densities and conditional expectations as 
well as its application to credit loss distributions with a partly continuous distribution. 

Section U introduces the model studied here and explains how to combine the kernel estimators 
from Section [3] with importance sampling for credit risk. 

Application of the algorithms introduced is illustrated with a numerical example in Section [5] 
We conclude with an assessment in Section [5] of what has been reached. 



2 Capital allocation 

In the following, we consider the following stochastic credit portfolio loss model: 

n 

?:=i 

Li, . . . ,L n > are random variables that represent the losses that a financial institution suffers on its 
exposures to borrowers i = 1, . . . ,n within a fixed time-period, e.g. one year. The random variable L 
then expresses the portfolio- wide loss. We denote by P[. . .] the real- world probability distribution that 
underlies model (|2.ip . In other words, P[. . .] is calibrated in such a way that it reflects as close as possible 
observed loss frequencies. P[L < £}, for instance, stands for the probability of observing portfolio- wide 
losses that do not exceed the amount I. The operator Ep[. . .] is defined as mathematical expectation 
with respect to probability P. In particular, Ep[X] reflects the real- world probability weighted mean of 
the portfolio-wide loss. 

As mentioned in the introduction, it is common practice for financial institutions to measure the risk 
inherent in their portfolios in terms of economic capital (EC). As credit risk, for most institutions, is 
considered to be most important, this is in particular relevant for credit portfolios. EC is commonly 
understood as a capital buffer intended to cover the losses of the lending financial institution with a high 
probability. This interpretation makes appear very natural the definition 

EC = VaR P . Q (L) -E P [L], (2.2a) 

where the Value-at-Risk (VaR) is given as a high-level (e.g. a = 99.9%) quantile of the portfolio-wide 
loss: 

VaR P , Q (i) = min{£ : P[L <£]> a}. (2.2b) 

2 Instead of trying to determine VaR contributions in a continuous or semi-continuous setting, some risk managers use 
contributions to expected short fall. This approach is very fruitful and has some other advantages (see [Kalkbrcnc r et all 
12004 iMerino and Nvfeleil 120041) . 
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Hence, if a financial institutions holds EC according to (|2.2a[) and charges the loans granted with upfront 
fees adding up to Ep[L], the probability that it will lose all its EC is not higher than 1 — a. Note that, 
despite its intuitiv e appeal, VaR as a risk m easure is criticised, e.g. for its potential lack of rewarding 
diversification (see lAcerbi and Tasche . 2002 . and the references therein). 



Active risk management involves more than just measuring portfolio- wide capital according to (|2.2ap . 
Additionally, it is of interest to identify which parts of the portfolio bind the largest portions of EC. 
The corresponding process of determining a risk-sensitive decomposition of EC is called capital alloca- 
tion. While for the expectation part Ep [L] of EC on the right-hand side of (|2.2a|) there is the natural 
decomposition 



EpLLl 



(2.3a) 



there is no such obvious decomposition 



VaR P , Q (L) = 5^VaR P , a (Z i |L) 



(2.3b) 



for the VaR-part of EC into risk contribution^. Indeed, the choice of the decomposition method depends 
on the conce pt of risk sens itivity adopted. Interpreting risk sensitivity as compatibility with portfolio 
optimization, iTasche (Il999h proved that the risk contributions VaRp. Q (Li | L) on the right-hand side of 
(I2.3b|) should be defined as directional derivatives, i.e. 



YsRp, a {Li | L) = 



dVaR PjQ (i + hLi) , 

dh ' /1=0 ' 



(2.4) 



As VaR is a positively homogeneoutQ risk measure, by Euler's theorem, then (|2.3b[) holds. (|2.4p displays 
the concept of risk contribution applied in this paper. Note however that in general, and in particular 
if th e distribution of L has no density, the derivative on the right-hand side (|2.4p need not exist. See, 
e.g., iTaschd ( 19991 Assumption (S)) for conditions ensuring existence of the derivative. Depending on 
the objective of the po rtfolio analysis , other approaches to determining risk contributions are reasonable, 
see, e.g., Section 3.1 of lTaschd (|2006h for an account of these. 



In general, no closed-form representations of VaRp ;Q (L) and the risk contributions VaRp :Q (L.; \L) are 
available. Therefore, often, these quantities can only be inferred from Monte-Carlo samples. This means 
essentially to generate a sample 



(L» ...,L«), t=l, 



T 



(2.5) 



and then to estimate the quantities under consideration on the basis of this sample. How to do this is 
quite obvious for VaR, but is much less clear for the risk contributions VaRp jQ ,(ii | L) as, in general, 
estimating derivatives of stochastic quantities without closed-form representation is a subtle issue. 

Fortunately, it turns out ( Gourieroux et al. . 2000l ; Lemus . 1999t Tasche . Il999h that, under fairly general 
conditions on the joint distribution of L and Li, the derivative (|2.4[) coincides with an expectation of the 
loss related to borrower i conditional on the event of observing a portfolio-wide loss equal to VaR. 



dVaR PjQ (£ + hLi) 
dh 



/i=0 = E P [L,|L = VaR P , Q (L)] 



(2.6) 



If P[L = VaRp iQ ,(L)] is positive, the conditional expectation on the right-hand side of (|2.6[) is given by 

Ep [Li l{L=VaR P , Q (L)}] 



E P [Li|L = VaR P , Q (L)] 



P[L = VaR PjQ (L)] 



(2.7) 



^Kalkbreneil d2005tl considers relation (12 . 3bl ) in a more general context. He calls it "linear aggregation". 
4 I.e. VaRp ta (hL) = hVaRp a (L) for positive h. 
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Even if P|X =_VaRp jQ (L)] is positive, its magnitude will usually be very small, such as 1 — a or less. 



Glassermanl (|2005l) shows how to apply importance sampling in such a situation in order to efficiently 
estimate EpfL^ | L = VaRp, Q (L)]. 



However, a crucial condition for (|2.6[) to hold exactly is the existence of a density of the distribution of 
L. The probability P[L = VaRp Q (L)] then equals zero, and consequently the right-hand side of (|2.7p is 
undefined^]. In thi s situation, th e conditional expectation Ep[L^ | L = VaRp , a (L)] is still well-defined (see, 
e.g., Remark 5.4 o f lTaschelll999h. but its e s timat ion from a sample like (|2.5[) requires more elaborated non- 
parametric methods. iMausser and Rosenl (|20Q4f ) suggest using an e stima tion method based on weighted 
combinations of order-statistics. We follow here Gourieroux et al. (2000) who applied kernel estimation 
methods for VaR contributions when optimizing returns in a portfolio of stocks. The kernel estimation 
procedures, however, have to be adapted to the r are-event characte r of credit risk. Therefore, in the 
remainder of the paper we modify the approach by iGourieroux et aT in a way that can be described as 
a combination of kernel estimation and importance sampling. 



3 Kernel estimators 

In this section, we introduce the classical Rosenblatt-Parzen kernel estimator for densities and the 
Nadaraya- Watson kernel estimator for conditional expectations in a way tha t links naturally t o the 
risk contribution concept of Section [2] The general reference for this section is Pagan and Ullah ( 19991 . 
Chapters 2 and 3). 

3.1 The Rosenblatt-Parzen kernel estimator for densities 

Assume that x%, . . . ,Xt is a sample of independent realisations of a random variable X with density /. 
The Rosenblatt-Parzen estimator fh with bandwidth h > for / can be constructed as follows: 

• Let X* be a random variable whose distribution is given by the empirical distribution corresponding 
to the sample x\, . . -,x T , i.e. P[X* = x t ] = 1/T, t=l,...,T. 

• Let £ a random variable with density (kernel) p. 

• Assume that X* and £ are independent. 

• Then the estimator fh is defined as the density of X* + h£: 

fh(x) = f h , Xu ..., XT (x) = 7TtEL^^)- (3 ' 1} 

If / and ip are appropriately "smooth" (see Pagan and Ullahl . 19991 Theorem 2.5 for details), it can be 
shown fc 
of /, i.e. 



T 1 — »oo T ->oo 

shown for h = Kt — — — > 0, HtT — — — > oo that fj lT (x) is a pointwisc mean-squared consistent estimator 



lim E[(/(x) - f hTiXl X T {x)Y] =0, (3.2) 

T — >oo 

with independent copies X\, . . . , Xt of X. While the Rosenblatt-Parzen density estimator is rather robust 
with respect to the choice of the kernel tp, it is quite sensitive to the choice of the bandwidth h. For the 
univariate case we consider here, efficient techniques like cross validation for the choice of the bandwidth 
are available. However, such more elaborated techniques usually involve some optimisation procedures 
that can be very time-consuming for large samples. As a consequence, for the purpose of this paper we 



5 This problem can be avoided by using the risk measure Expected Shortfall (sec, e.g., Accrbi and Taschc, 2002) instead 
of VaR. With the definition of Expected Shortfall slightly simplified for practical purposes, 112.71 1 then reads Ep[L; \L > 
VaR P , a (L)] = (l-Q)- 1 Ep[L 1 l { L>v a Rp. Q (L)}]- 
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confi ne ourselves to applying a simple rule of thumb by Silverman (cf. Chapter 2 in iPagan and Ullahl . 
19991) 

h = 1.06cr T~ 1/5 , (3.3) 

where a denotes the standard deviation of the sample x\, . . . ,Xt- Moreover, we choose the standard 
normal density as the kernel ip. 



3.2 The Nadaraya- Watson kernel estimator for conditional expectations 



Assume that (x\, yi), . . . , (xt, Ut) is a sample of realisations of a random vector (X, Y) where X has a 
density /. The Nadaraya- Watson estimator E^[F|X = x\ with bandwidth h for E[Y | X = x] can be 
constructed as follows: 



• Let {X* , Y*) a random vector whose distribution is given by the empirical distribution correspond- 
ing to the sample {xi,y{), (x T ,yr), i.e. P[(X*,Y*) = (x t ,y t )} = 1/T. 

• Let £ a random variable with density (kernel) p>. 

• Assume that (X*,Y*) and £ are independent. 

• Then the estimator Eh[Y | X = x] is defined as the expectation of Y* conditional on X* + h£ = x: 

E h [Y\X = x] = V hAxi , Vl) _ (XTtVT) [Y\X = x] = ^T Vt ^j£) ■ ( 3 - 4 ) 

If / and p are appropriately "smooth" (see Pagan and Ullah . 19991 Theorem 3.4 for details), it can be 
shown for h = hr T ^ co > 0, hrT T ^°° > 00) anc i f(x) > that E/j[Y | X = x] is a pointwise consistent 
estimator of E[Y | X = x], i.e. 

lim P[ | E[Y | X = x] - % lT , ( x 1 ,x 1 )....,(x T ,x T )[ Y \X = x]\>e] =Q, e > arbitrary, (3.5) 

with independent copies (X\, Yi), . . . , (Xt, Yp) of (X,Y). The construction of the Nadaraya- Watson 
estimator (|3.4p as described above allows to interpret the conditional expectation estimation problem as 
an extended density estimation problem. This suggests to choose the same bandwidth h and the same 
kernel <p for the estimators (|3.ip and (|3.4p . 



Remark 3.1 Assume that in the random vector (X,Y) the X -component is a sum of random variables 
X\, . . . , X n and that we are interested in estimating E[X, \X = x], i = 1, . . . , n. Define (X*, . . . , X*), 
analogously to (X*,Y*) as the "empirical" version of (Xi, . . . , X n ). According to <\'SA\i , the Nadaraya- 
Watson estimator o/E[Xj | X = x] can then be specified as 

E h [Xi\X = x] = E[X;\X* + ht = x], (3.6a) 

with an appropriate auxiliary variable £ independent of (X^, . . . , X*). If the same bandwidth h is applied 
for all i, then from representation (|3.6a[) follows 

n 

Y J V h [X l \X = x] = x-hE[£\X* + h£ = x]. (3.6b) 

i=l 

As we will note in Section \ 5.S\ when commenting on Tabled the size of the difference of the left-hand 
side of (|3.6b[) and x can be regarded as providing additional information on the choice of the bandwidth 
h. We will apply the multiplicative adjustment suggested bv lEvverlein and Smillie\ kOOd , Equation (7)) 
to force additivity in the sense of (|2.3b[) on the estimated contributions to VaR. 
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3.3 Application to credit losses 



For every credit risk portfolio, there is some positive probability of observing no losses. While for large 
portfolios, this probability will usually be negligibly small, we will see by the example from Section [5] 
that the probability of zero loss can be of significant magnitude for smaller portfolios. As a consequence, 
we cannot assume that the loss variable L from model (|2.1[) has an unconditional density. For otherwise 
P[L = 0] would be zero. Hence, at first glance, applying the estimators (|3.1[) and (|3 .4[) seems not possible 
in the context of model (|2.ip . The following assumption, however, allows us to deal with this problem. 



Assumption 3.2 There is a random vector (Si, . . . ,Sk), (called systematic factors,), with the following 
properties: 



(i) The loss variables Li in (|2.1|) are independent conditional on realisations of (Si, . . . , Sk). 

(ii) For each i = 1, . . . ,n, there are probabilities Pi(si, . . . , Sk) G [0, 1] and densities < fi(l, si, . . . , Sk) 
such that the distribution function of Li conditional on (Si, . . . , Sk) is given for £ > by 

P[L i <e\(Si,...,Sk) = (s 1 ,...,s k )} = 

l-Pi(si,...,Sk)+Pi(si,...,Sk) fi(x,si,...,Sk)dx. (3.7) 



Under Assumption 13. 21 it is easy to derive the following result on the representation of the unconditional 
distribution of the portfolio-wide loss L. 

Theorem 3.3 Define I n = {l,...,n} and write (^) for the multiple convolution of densities. Then, 
under assumvtion \3.Sl for £ > the distribution function of the loss variable L from (|2.1|) can be written 
as 

P[L<e]=p+(l-p) [ f(x)dx, (3.8) 
Jo 

with 

p = P[Li=0,...,L n = 0], 
f(x)= ]T E P [H Pi (Si,...,S k ) J] (l-p i (S 1) ...,S k ))((g)f i (.,S 1 ,...,S k ))(x) . 



Although, due to the involved multiple convolutions, (|3.8p is not really useful for calculating the distribu- 
tion of L, it allows us to assume that, conditional on being positive, the portfolio-wide loss has a density, 
i.e. for £ > 



P[L< £ \L > 0] = / f(x)dx. 

Jo 



(3.9) 



Define P* by 

P*[A] = P[A\L > 0] (3.10) 
for any relevant event A. The following lemma is then obvious. 



Lemma 3.4 If the probability P* is given by (|3.10[) . then for £ > the expectations conditional on L 
with respect to P* and P are identical. In particular, in the context of model (|2.1[) for all i = 1, . . 
and £ > we have 

E P [L l \L = e] = E P .[Li\L = £\. 
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Note that a sample (|2.5p generated under measure P becomes a sample generated under measure P* 
when all (n + l)-tuples (L^\ . . . , L$ ) with = are eliminated. For this sub-sample then 
the preconditions for applying estimators (|3.ip and (|3 .4[) are satisfied. This observation leads to the 
following algorithm for estimating the risk contributions according to (|2.6p for model (|2.1[) by Monte 
Carlo sampling. 

Algorithm 3.5 

i. Generate a sample like (|2.5p /rom £/ie real-world probability measure P. 

Determine an estimat^ £ of VaRp jQ ,(L) /rom i/iis sample. 
3. Extract the sub-sample with L^' > /rom the previous sample. 

4- Calculate, on the basis of the sub-sample with > 0, the bandwidth h* according to (|3.3D . 

5. Calculate the estimates for Ep[Li \ L = VaRp jQ ,(i)], i = 1, . ..,71, on the basis of the sub-sample 
with Z,W > 0, according to (|3T3|) as E ft . [£, | i = I] (cf. Remark \3~f\) . 

In the following section, we will modify this algorithm by incorporating importance sampling for reducing 
the variances of the estimates. 



4 Importance sampling for credit risk 



McNeil et ai1(|2005l ), at the beginning of Chapter 8.5, write "A possible method for calculating risk mea- 



sures and related quantities such as capital allocations is to use Monte Carlo (MC) simulation, although 
the problem of rare event simulation arises." They then explain, in the context of risk contributions to 
Expected Shortfall, that "the standard MC estimator . . . will be unstable and subject to high variabil- 
ity, unless the number of simulations is very large. The problem is of course that most simulations are 
'wasted', in that they lead to a value of L which is smaller than VaRp jQ ,(L)." What applies to contribu- 
tions to expected shortfall applies even more to VaR contributions as these, according to (|2.6p . are related 
to events still more rare and actually in our case of probability zero. Thus, the key idea with importance 
sampling is to replace the real-world probability measure P by a probability measure Q which puts more 
mass on the interesting events. 

Some technical assumptions are needed to guarantee that such a replacement of probabilities does really 
work. 

Assumption 4.1 

• The probability measures P and Q are defined on the same measurable space (fi,^ 7 ). 

• There is a measure (i on (f2, J 7 ) such that both P and Q are absolutely continuous with respect to 
\i. Denote by f the density o/P and by g the density o/Q. 

• / > implies g > 0, i.e. P is absolutely continuous with respect to Q. 



Under Assumption 14. 11 the likelihood ratio 

R = - (4.1) 
9 

is Q-almost surely well-defined. This implies that any expectation with respect to P can be expressed as 
an expectation with respect to Q, i.e. for any integrable X holds 

EpLY] = E Q [RX]. (4.2) 



"We may assume £ > as for real- world portfolios the case I = seems very unlikely. 
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However, according to (|2.6p . for the purpose of this paper it is a conditional rather than an uncondi- 
tional expectation we are interested in. The following proposition gives a result analogous to (|4.2p . for 
conditional expectations. 



Proposition 4.2 Let P and Q be probability measures as in A s sumption \4 ■ 1\ with [i-densities f and g 
respectively. Define the likelihood ratio R by (|4.ip . Then we have for any sub-o -algebra A of J 7 and any 
integrable random variable X 



Proof. See e.g. iKlebanerl (|2005l . Theorem 10 



□ 



The choice of an appropriate measure Q for use in (|4.2|) and (|4.3j) is not at all obvious. Inspecting (|2.1| 
and (|2.6p . it becomes clear that, for the purpose of this paper, we are interested in events that result in 
losses clos(0 to 1 = VaRp. Q (L), as defined by (|2.2bp . 

In Section[5l we will specifically consider a numerical example that satisfies the following conditions which 
reflect assumptions commonly made for models in industry. 



Assumption 4.3 The portfolio-wide loss L is given by ()2.ip . L\,... ,L n are given as 

U = Ail Di . 



D x ,. 
P[D, 

of the Di 



D n are independent (default) events, conditional on a set of systematic factors (Si, . . . , Sfc), with 
Pi . The loss severity variables A±, . . . , A n are positive and independent, as well as independent 
, Sfc. 



, D n and the Si , . 



The distribution of Ai is specified via its density a,(s) > 0. 



Note that, under the condition of the Ai having densities, Assumption 14. 31 is a special case of Assumption 
3.21 As a consequence of this, according to Theorem 13.31 the loss variable L from (|2.ip has, in theory, a 
density for its positive realisations. 

Merino and Nvfeler ( 2004 ) and Glasserman and Ll ( 2005 ) suggest a nested simulation procedure where 
exponential twisting is applied for the estimation of expectations conditional on the systematic factors. The 
resulting conditional loss distribution then can essentially be described again by Assumption 14. 3[ with 
independence instead of conditional independence and modified probabilities of default. Unfortunately, 
this nested procedure cannot be applied for our problem of estimating Ep [Li\L = ■] as there is no 
independence conditional on L. In general, we have a(L) <jt a(Si, . . . , St) and a(Si, . . . , Sfe) <£. a{L). 
Therefore, no nesting of conditioning is applicable either. Exponential twisting can be applied nevertheless 
but does not yield satisfactory results as it is not clear how an optimal tilting parameter should be 
determined. 



Approaches by iKalkbrener et alj (|2004f ) and iGlasserman and Lil (120051 ) are promising alternatives to 
exponential twisting. In these approaches, the importanc e sampling measu r e Q i s created by changing 
the means of the systematic factors from Assumption 14.31 IKalkbrener et all (|2004l ) suggest to determine 
the factor means appropriate for importance sampling by s olving a minimisation problem for the variance 
of the estimator they consider. IGlasserman and M (120051) instead look for factor means that make the 
mode of the factor distribution coincide (approximately) with the mode of the "zero- variance importance 
sampling" distribution. However, these approaches are more complex compared to exponential twisting 
as they involve - in the case of a multi-factor model - choosing several parameters instead of only one as 
required by exponential twisting. 



The approach we follow in this paper is close in spirit to the approach by IGlasserman and Li ( 2005 ). 
Define the importance sampling probability measure Q = Q M by Assumption 14.31 but replace the vector 



7 In general, I itself will have to be estimated. Thus, at first glance, it seems strange to choose it as a basis for finding 
Q. However, in a first step, for instance, it can be replaced by a rough estimate and be refined in further stages of the 
estimation procedure. 
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(Si, ■ ■ ■ ,Sk) of systematic factors by a vector (Sf, with 

5; = Si -E P [$] (4.4a) 

where fi = (fii, . . . , fj>k) satisfies 

Hi « E P [5 i |L = VaR P>a (L)]. (4.4b) 
The conditional expectations in (|4.4b|) are estimated by applying the Nadaraya- Watson estimator (|3.4j) . 



Algorithm 13.51 for estimating the risk contributions according to (|2.6[) for model (|2.1|) by Monte Carlo 
sampling has to be modified as follows when importance sampling is applied. 



Algorithm 4.4 

1. Generate a sample (L^\ , . . . , L„ ), i = 1, . . . , Ti, from the original sampling probability measure 
P. 

2. Determine an estimate £ of VaRp iQ (T) /rom this sample. 

3. Calculate, on the basis of the sub-sample with > 0, i/ie standard MC simulation bandwidth h* 
according to (|3.3j) . 

^. Estimate the mean shift parameters Hi according to (|4.4b[) /rom the sample (Lft\ L\ , . . . , T„ ), 
i = 1, . . . , Ti, /or instance by applying the Nadaraya- Watson estimator. 

5. Generate a sampl^ (L^\ Ti , • ■ ■ , T^, iiW), t = 1, . . . , T2, /rom f/ie importance sampling proba- 
bility measure Q M . By' denotes realisations of the likelihood ratio as defined by (|4.1D . 

6. (Optional) Determine a refined estimat^\ i o/VaRp. Q (T) from this sample. 

7. Calculate, on the basis of the sub-sample with LW > 0, the importance sampling bandwidth h* 
according to (|3.3p . 

8. Estimate, on the basis of the sub-sample with > 0, the Q M - conditional expectations Eq [i2 | T = 
£] and Eq (j [i?T i | T = £], i = 1, . . . , n according to p.4p wii/i bandwidth h* . 

9. Calculate the estimates for Ep[T^ | T = VaRp. Q (T)], i = 1, . . . ,n according to (|4.3[) . inserting the 
estimates from the previous step. 



5 Numerical example 



This section is divided into two parts. In Sub-section 15.11 we provide details of the portfolio model that 
is used for the simulation study and describe the numerical results we seek to obtain. In Sub-section 15.21 
we comment on the results and point out their essential features. 



5.1 Description of simulation study 



We examine the performance of importance sampling estimation of VaR contributions under Assumption 
14.31 Due to restriction in computational power, the portfolio we consider is relatively small, with 96 assets. 
Note that, with such a portfolio size, there may be a significant positive probability of not observing any 
losses, demonstrating that the issue tackled in Lemma 13.41 has some relevance. 

8 The new sample need not necessarily be generated by a new Monte Carlo simulation. Alternatively, as done in Section 
[5] the new sample can be created from the previous sample by substituting (SJ, . . . , S 1 ^) from ]4,4aH for (Si, . . . , S^). In 
this case T2 = Ti. 

9 For instance, by ordering the pairs (I", JiW) in descending order according to the L-component and selecting the 
largest t with R{k) 

< 1 — a. Take then ZA*' as the estimator. 
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The dependence structure of the portfolio is determined by four correlated systematic factors (Si, S2, S3, 5*4 ) 
which are each standard normal and are jointly normally distributed with correlation matrix 



/ 1 


0.75 


0.05 


0.75 


1 


0.05 


0.05 


0.05 


1 


\0.05 


0.05 


0.25 



0.05\ 

0.05 

0.25 

1 1 



(5.1) 



According to (|5.1|) . factors Si and S2 are strongly correlated, factors S3 and S4 are moderately correlated, 
and the pairs (Si, S2) and (S3, S4) are weakly dependent. 

Each factor corresponds to one sector that includes 24 assets. In each of the four sectors, the risk 
characteristics of the 24 assets are specified as shown in Table [TJ According to Table [TJ each of the 
four sectors in the portfolio includes 12 high (2%) PD (probability of default) assets and 12 low (0.5%) 
PD assets. In both of these two sub-sectors there are 4 high exposure (25$), 4 medium-size exposure 
(5$), and 4 low exposure (1$) assets. Of the 4 high exposure assets, all have equal LGD (loss given 
default) mean 50% but 2 assets have high LGD variance (12.5%) and 2 assets have low LGD variance 
(3.125%). Similarly, among the 4 medium-size (low) exposure assets of equal LGD mean, there are 2 high 
LGD variance and 2 low LGD variance assets. Note that for each combination of sector, PD, exposure, 
and LGD variance there are two assets with identical risk characteristics. This feature of the portfolio 
composition is intended to deliver a rough assessment of estimation uncertainty due to the Monte Carlo 
simulation. For, by symmetry, risk contributions for assets with identical risk characteristics should be 
equal but will not be when estimated by Monte Carlo simulation. 

Assumption 5.1 For any asset i £ {1, 2, . . . ,96} in the portfolio, its individual loss variable Li in the 
sense of Assumption \4~^\ is specified by the asset's sector k(i) € {1,2,3,4}, probability of default PDi, 
exposure Vi, mean loss given default LGDi, and LGD variance varLGDi (as shown in Table[lty as follows: 



The default event of the asset is given by D L = S^^i) +y/l — r£i < <5>~ 1 (PDi)} where <f> denotes 
the standard normal distribution function, r = 0.18 (equal for all i), and ^1,^2, • ■ ■ ,£96 « r e i-i.d. 
standard normal. 

The loss severity variable is given by Ai = Vi Bi where By , B% , . . . , Bqq are independent beta- 
distributed random variables with E[£?i] = LGDi and var[i?i] = varLGDi. 



The constant r in the definition of the default events is the loading of the systematic risk in this model. 
Its value was chosen with a view on th e asset correl ations in the Basel II corporate risk weight formula 
which have a range from 0.12 to 0.24 (jBCBSl . I2OO6I paragraph 272). The loss severity distributions of 



the high LGD variance assets in Assumption 15.11 are U-shaped (beta parameters a = 0.5 and b = 0.5), 
the severity distributions of the low LGD variance assets are bell-shaped (beta parameters a = 3.5 and 
b = 3.5). 

By portfolio construction, on the one hand Sectors 1 and 2 are identical with respect to their composition 
and risk characteristics. On the other hand, this holds also for Sectors 3 and 4. As Sectors 1 and 2 are, 
however, stronger correlated than Sectors 3 and 4, one should expect that due to this concentration the 
risk contributions of assets in Sectors 1 or 2 arc higher than the risk contributions of the corresponding 
assets in Sectors 3 or 4. 

The simulation excrcisd^l we conduct is structured as follows: 

1. We run 25 times a Monte Carlo simulation with 50,000 joint realisations of the losses of all the 
assets in the portfolio. 

2. In each simulation run, the following quantities are estimated: 



'The R-script for the calculations can be down-loaded at http://www-m4.ma.tum.de/pers/tasche/. 
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• Portfolio-wide VaR at 99.9% level (by standard MC and importance sampling MC, sec Table 
[5] for the results). 

• For each factor, the conditional expectation Ep[5j | L = VaRp :Q (L)] of the factor conditional 
on L = VaR P , Q (£) (by standard MC, Table©. 

• For each sector, sector-stand-alone VaR at 99.9% level (by standard MC and importance 
sampling MC, Table EJ). 

• For each asset, the contribution of the asset to portfolio-wide VaR at 99.9% level (by standard 
MC and importance sampling MC, Table [7]). Additionally, for each asset, the contribution of 
the asset to portfolio- wide VaR at 99.9% level by importance sampling MC with reduced and 
enlarged respectively bandwidth for the conditional expectation, Table [5]) . 

3. Based on the estimates of step 2, in each simulation run the following quantities are calculated: 

• The portfolio-wide mean of the loss, standard deviation of the loss, and the probability of 
observing a los^3 (Table 5|) . 

• The ratio of the sum of the VaR contributions of all assets and portfolio- wide VaR (Table [6]). 

• For each sector, the contribution of the sector to portfolio-wide VaR at 99.9% level (as the 
sum of the VaR contributions of the assets in the sector, Table©. 

• The portfolio- wide diversification indeJ^I (Table [T0|) 

KL) = V^W-EW (5.2a) 
Er=i(VaR P , a (ii)-E[L l ]) V j 

• For each asset i, the marginal diversification index (Table [TT]) 

UT In VBR Pia (L i \L)--E[L i ] 

l(LtlL) = VaR P ,(L i )-E[L i ] ' ^ 

For each sector, its marginal diversification index (defined as the sum of the VaR contributions 
of the assets in the sector divided by the VaR of the sector, Table [T2|) . 



Marginal diversification indices as defined by (|5.2bj) represent a direct application of risk contri- 
butions for risk concentration analysis. By construction, I(Li | L) > I(L) imp lies that reduction 
of the exposure to asset i will improve portfolio diversification (|Taschd . I2006L Section 4). Here, 



"diversification" is understood in a relative sense, namely comparing the actual economic capi- 
tal assigned to the portfolio under consideration to the economic capital assigned to a worst case 
portfolio composed of co-monotonic loss variables. 

4. We report three results for each of the values that is estimated or calculated from estimates: 

• The result of the first simulation run. 

• The mean of the results of all 25 simulation runs (as approximation of the true value). 



The coefficients of variational of the results of all 25 simulation runs (as measure of estimation 
uncertainty) . 



11 For further reduction of the estimation variance, the standard deviation and sample size of the positive losses for 
the bandwidth (both for standard MC as well as for importance sampling) according to 1 13. 31 1 are not estimated but 
numerically calculated. Under Assumption 15.11 this can be done exactly for the standard deviation of the positive losses 
and approximately for the sample size of the positive losses, by approximating the distribution of the number of defaults 
via mo ment ma t ching by a negative binomial distribution. 

12 Scc Taschc (2006, Section 4) for a motivation of this definition and the definition of marginal diversification indices. 
The diversification indices are calculated on an unexpected loss (UL) basis, in accordance with definition 112. 2al of economic 
capital. Note that the value of the portfolio-wide diversification index depends upon whether the portfolio is decomposed 
into assets or into sectors since in the latter case the diversification potential is larger. 

13 The coefficient of variation of a sample is defined as the ratio of the sample standard deviation and the sample mean. 
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5.2 Comments on the results 



Table [2j With respect to the estimation of portfolio-wide VaR, according to the results displayed in 
Table O there is a minor ad vantage in using importanc e sampling. With multi-step importance sampling 
as suggested for instance by iGlasserman and Lil (|2005l ) this advantage could be increased. However, the 
purpose of this paper is to deal with the estimation of VaR contributions. Therefore, here we need not 
look into the details of efficient estimation of VaR itself. As indicated by the relatively small sample 
variation coefficients both for standard as well as importance Monte Carlo sampling, the estimates in the 
first simulation run do not differ much from the means of 25 simulation runs. 



Table \3[ Table [3] shows that the estimates for the shift constants of the distribution of the systematic 
factors according to (|4.4aj) and (|4.4b|) arc not very stable. This is indicated both by the high sample 
variation coefficients as well as by the 1st run estimate of the constant for the fourth factor which differs 
much from the sample mean. As the strongly correlated factors 1 and 2 might contribute more to portfolio 
risk as measured by VaR than the factors 3 and 4, it is no surprise that the indicated shifts of factors 1 
and 2 are larger than the ones of factors 3 and 4. 



Table [4j The difference between the loss distribution under the original probability measure P and 
the importance sampling measures is demonstrated by Table [U The mean loss grows by more than ten 
times and then slightly overshoots the portfolio VaR under the original probability measure P. The loss 
standard deviation increases by four times. The proportion of portfolio loss realisations with positive 
losses grows from 60% to almost 100%. Hence the proportion of the sample that can be used for kernel 
estimation is much greater in the case of importance sampling. Note that the characteristics of the shift 
loss distribution do not seem to vary much in the 25 simulation runs. 



Table [5j By construction, from a risk perspective the four portfolio sectors are identical when considered 
stand-alone. This is confirmed by the estimates of stand-alone sector VaR as displayed in Table O The 
significantly higher sample variation coefficients of sectors 3 and 4 in the case of importance sampling 
show that the shift factor distribution according to (|4.4ap and (|4.4b[) is not very well suited for the 
estimation of stand-alone sector characteristics. 



Table [6j As indicated in Remark 13.11 although in theory the sum of the VaR contributions according 
to (|2.6|) equals VaR, a sum of VaR contribution estimates made by kernel estimation can differ from VaR. 
Table |6] displays for different estimation approaches how large the difference can be. In general, it is larger 
for importance sampling and increases with the kernel estimation bandwidth. As we will see, the choice of 
the importance sampling bandwidth according to Silverman's rule of thumb (|3.3p seems to be a reasonable 
compromise between reduction of sample variation by oversmoothing and unbiascdness as measured by 
the difference between the sum of VaR contributions and VaR. Note that, to make comparable the VaR 
contribution estimates by different approaches, the contributions as displayed in Tables \7\ [SJ and O are 
normalised such that their sum equals portfolio VaR. 



Tables [7] and [SJ As by symmetry results for sectors 2 and 4 are not essentially different from the 
results for sectors 1 and 3 respectively, Tables [7] and [8] display the results of VaR contribution estimates 
at asset level only for sectors 1 and 3. The tables allow to compare the estimation performance as 
yielded by four different approaches: kernel estimation based on standard Monte Carlo sampling, kernel 
estimation based on importance sampling with bandwidth chosen according to Silverman's rule of thumb, 
kernel estimation based on importance sampling with reduced bandwidth, and kernel estimation based 
on importance sampling with enlarged bandwidth. The average results of all the approaches do not differ 
too much - which is also a consequence of the normalising applied to have the sum of the contributions 
equal to portfolio VaR. The variation of the estimates, as indicated by their coefficients of variation, 
is clearly highest for standard Monte Carlo, followed by importance sampling with 50% of Silverman's 
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bandwidth. Importance sampling with 200% of Silverman's bandwidth has lower sample variation than 
importance sampling with 100% of Silverman's bandwidth. As noticed above, the latter, however, need 
not be so much adjusted as the former. For this reason, importance sampling with bandwidth according 
to Silverman's rule might be considered the best estimation method for VaR contributions among the 
four approaches discussed here. 

Taking into account that the values in Tables [7] and [5] were calculated on the basis of 25 simulation runs 
each with 50,000 loss realisations, the results are somewhat disappointing as the sample variation is still 
large in particular for small exposures with low probability of default. As demonstrated by Figure [TJ 
there is nevertheless a clear gain in estimation efficiency by the application of importance sampling. 
Additionally, in contrast to the standard Monte Carlo approach, the importance sampling approaches 
have no problem with seeming zero VaR contributions of some assets (compare the first run results) . 

Note that the correlation structure of the portfolio is clearly reflected in the VaR contributions as the 
contributions by assets in sector 1 are significantly higher than the contributions by the assets with 
similar risk characteristics in sector 3. Note also that LGD variance has a strong impact on an asset's 
VaR contribution, in particular for those assets with the high exposures. The greater the LGD variance, 
the greater the VaR contribution. 

Table [9J Tabic [5] displays at sector level what Table [7J shows at asset level. Compared to the asset 
level, sample variation at sector level is slightly less. Again, by importance sampling there is some gain in 
estimation efficiency compared to standard sampling. Note in particular for sector 1 the misleading first 
run standard estimate of the VaR contribution which seems to indicate that sector 1 is less correlated to 
the rest of the portfolio than sectors 2, 3, or 4. 

Table 1101 As demonstrated in Tabic [101 estimates of portfolio- wide diversification indices in asset and 
sector context are fairly stable, and there is not much difference in efficiency between standard and 
importance sampling. The specification of the context in which the diversification indices are calculated 
indicates the scope of possible actions for a reduction of portfolio concentration. Sector context means 
that only the relative weights of the sectors may be changed but not relative weights of single assets. 
Superficially, the sector diversification index looks worse. This is caused by the fact that its denominator 
- compared to the asset context - is less because the sector VaR figures already incorporate a lot of 
diversification. In general, it does not make sense to compare diversification indices that were calculated 
in different contexts. Portfolio- wide diversification indices should rather be compared to the corresponding 
asset or sector diversification indices because this way guidance can be provided on how to change the 
portfolio for better diversification. 

Table 1 1 11 Not surprisingly, as Tabic [Tl] shows importance sampling estimates are more efficient than 
standard sampling estimates also for the estimation of asset-level marginal diversification indices. In 
particular, importance sampling avoids observing negative diversification indices even in the first run 
estimates where estimation is much more uncertain than in the "mean of all runs" columns. Note 
however, that also the importance sampling first run estimates arc misleading in so far as they seem to 
indicate that asset 4 (with low LGD variance) in sector 1 is the most dangerous in the portfolio. In fact, 
assets 1 and 2 in this sector with same exposure size, PD, and mean LGD are more dangerous - as is 
correctly shown in columns 4 and 7 of Table ITTI - because they have got higher LGD variance. Note also 
that exposure concentration can "override" concentration caused by correlation. This is illustrated by 
assets 1 to 4 of sector 3 whose diversification indices are higher than the portfolio-wide index although 
sector 3 clearly has less correlation to the rest of the portfolio than sectors 1 or 2. 

Table 1121 Also at sector level the importance sampling estimates of the marginal diversification indices 
display less variation than the standard sampling estimates. In particular, the standard sampling first run 
estimate of the index for sector 1 is lower than the portfolio-wide diversification index. The misleading 
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conclusion could be that shifting weight to sector 1 contributes to portfolio diversification, in spite of the 
strong correlation between sectors 1 and 2. In contrast, the importance sampling first run estimate of the 
diversification index for sector 1 correctly indicates that the index is larger than the portfolio- wide index. 
However, it also indicates erroneously that sector 1 is worse for portfolio diversification than sector 2. 



6 Conclusions 

In general, determining VaR contributions in a credit portfolio risk model that involves continuous loss 
given default rate distributions is a non-trivial task. In the context of the common approach by means 
of Monte-Carlo simulation, we have discussed how to adapt kernel estimation methods for this problem 
and how to combine them with importance sampling. Importance sampling, in the form of a shift of the 
distribution of the systematic factors, is applied here since the variability of the estimates is quite strong, 
as a consequence of the rare-event character of credit risk realisations. 

The numerical example presented in Section [5] illustrates that the gain in estimation efficiency by these 
methods is significant. It also reveals, however, that the results yielded with these methods are not yet 
too satisfactory in so far as in particular the variability of estimates of VaR contributions for exposures 
with very small PDs re mains still quite large. It seems worthwhile to analyse in more detail how multi- 
step approaches e.g. by asserinan et all (|2005h have to be modified for successful application on the 



estimation of VaR contributions as studied here. Further research on this issue could be useful. 
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Table 1: Risk characteristics of assets. Identical for all four sectors. 



Asset No. 


PD 


Exposure 


LCD mean 


LCD variance 


Stand-alone VaRgg 9 % 


1 


0.02 


25 


0.5 


0.125 


24.846 


2 


0.02 


25 


0.5 


0.125 


24.846 


3 


0.02 


25 


0.5 


0.03125 


19.778 


4 


0.02 


25 


0.5 


0.03125 


19.778 


5 


0.02 


5 


0.5 


0.125 


4.969 


6 


0.02 


5 


0.5 


0.125 


4.969 


7 


0.02 


5 


0.5 


0.03125 


3.956 


8 


0.02 


5 


0.5 


0.03125 


3.956 


9 


0.02 


1 


0.5 


0.125 


0.994 


10 


0.02 


1 


0.5 


0.125 


0.994 


11 


0.02 


1 


0.5 


0.03125 


0.791 


12 


0.02 


1 


0.5 


0.03125 


0.791 


13 


0.005 


25 


0.5 


0.125 


22.613 


14 


0.005 


25 


0.5 


0.125 


22.613 


15 


0.005 


25 


0.5 


0.03125 


16.51 


16 


0.005 


25 


0.5 


0.03125 


16.51 


17 


0.005 


5 


0.5 


0.125 


4.523 


18 


0.005 


5 


0.5 


0.125 


4.523 


19 


0.005 


5 


0.5 


0.03125 


3.302 


20 


0.005 


5 


0.5 


0.03125 


3.302 


21 


0.005 


1 


0.5 


0.125 


0.905 


22 


0.005 


1 


0.5 


0.125 


0.905 


23 


0.005 


1 


0.5 


0.03125 


0.66 


24 


0.005 


1 


0.5 


0.03125 


0.66 
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Tabic 2: Standard Monte Carlo and importance sampling estimates of portfolio VaR at 99.9% level 
and coefficients of variation of estimates. 



Sampling method 


1st run 


Mean of all runs 


Coef. of var. 


Standard MC 


69.29 


69.22 


1.88% 


Imp. samp. MC 


68.4 


68.7 


1.22% 



Tabic 3: Standard Monte Carlo estimates of conditional expectations of systematic factors condi- 
tional on "Loss equals portfolio VaR at 99.9% level" and coefficients of variation of estimates. 



Factor 


1st run 


Mean of all runs 


Coef. of var. 


1 


-1.327 


-1.449 


16.38% 


2 


-1.313 


-1.471 


14.23% 


3 


-0.852 


-0.854 


22.04% 


4 


-1.244 


-0.87 


34.15% 



Tabic 4: Portfolio loss distribution characteristics mean of loss, standard deviation of loss, and 
probability of observing positive losses for original and importance sampling distribution. Charac- 
teristics for original distribution are calculated without simulation only once before the 1st simula- 
tion run. 



Sector 


Original distribution 


Importance sampling distribution 






1st run 


Mean of all runs 


Coef. of var. 


Expected loss 


6.2 


70.101 


73.204 


15.73% 


Stddcv of loss 


10.359 


41.406 


42.204 


8.14% 


Prob. pos. loss 


59.411% 


99.911% 


99.881% 


0.15% 



Tabic 5: Standard Monte Carlo and importance sampling estimates of stand-alone sector-VaR at 
99. 9%o level and coefficients of variation of estimates. 



Sector 


Standard MC 


Importance sampling 




1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Coef. of var. 


1 


39.434 


41.379 


2.79% 


40.366 


41.253 


1.95% 


2 


42.302 


41.524 


2.67% 


39.969 


41.487 


2.53% 


3 


40.954 


41.504 


2.49% 


40.369 


41.211 


3.37% 


4 


41.15 


41.285 


2.08% 


42.674 


41.623 


4.71% 



Tabic 6: Ratio of sum of asset VaR contributions and portfolio-wide VaR (at 99.9% level) for 
Standard Monte Carlo, importance sampling, importance sampling with 50% of bandwidth, and 
importance sampling with 200% of bandwidth. 



Sampling method 


1st run 


With mean contributions and VaRs of all runs 


Standard MC 


99.69% 


99.6% 


Imp. samp. MC 


96.84% 


96.2% 


50% bandwidth 


99.15% 


98.96% 


200% bandwidth 


85.31% 


85.29% 
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Table 7: Standard Monte Carlo and importance sampling estimates of asset contributions to portfolio-wide 
VaR at 99.9% level and coefficients of variation of estimates. Only sectors 1 and 3. See Figure^ for a 
graphical comparison of the coefficients of variation. 



Sector 


Asset 


Standard MC 


Importance sampling 






1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Coef. of var. 
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3 73991 


4.13472 


39 57% 
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Figure 1: Graphical comparison of coefficients of variation from Table^ 



Sector 1 : Coefficients of variation of VaR contribution estimates 
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Sector 3: Coefficients of variation of VaR contribution estimates 
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Table 8: Importance sampling estimates with half and double of bandwidth according to (|3.3[) for asset 
contributions to portfolio-wide VaR at 99. 9% level and coefficients of variation of estimates. Only sectors 
1 and 3. 



Sector 


Asset 


Imp. samp., 50% bandwidth 


Imp. samp., 200% bandwidth 






1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Coef. of var. 


1 


1 


3.26946 


4.29182 


38.02% 


4.05739 


4.00807 


16.07% 


1 


2 


4.66724 


3.94153 


16.02% 


4.03265 


4.09244 


12.27% 


1 


3 


1.98064 


2.64532 


19.71% 


2.14906 


2.72078 


16.23% 


1 


4 


3.80636 


2.57832 


23.22% 


3.34566 


2.8715 


18.5% 


1 


5 


0.26897 


0.26431 


41.28% 


0.28143 


0.25674 


16.95% 


1 


6 


0.33207 


0.29524 


34.96% 


0.30172 


0.26413 


12.17% 


1 


7 


0.26506 


0.23521 


19.71% 


0.24753 


0.2373 


16.24% 


1 


8 


0.24949 


0.24084 


25.59% 


0.23923 


0.24002 


14.17% 


1 


9 


0.05815 


0.04473 


34.66% 


0.04556 


0.04185 


15.54% 


1 


10 


0.05494 


0.03912 


21.92% 


0.04404 


0.04157 


10.76% 


1 


11 


0.04378 


0.04123 


20.37% 


0.04292 


0.04141 


12.44% 


1 


12 


0.04586 


0.04448 


38.31% 


0.04482 


0.04233 


12.57% 


1 


13 


1.34413 


1.3799 


43.47% 


1.17617 


1.26651 


18.33% 


1 


14 


1.90522 


1.31631 


31.05% 


1.56665 


1.21016 


17.01% 


1 


15 


0.84645 


1.04631 


45.4% 


0.85945 


0.88251 


20.24% 


1 


16 


1.27248 


0.87883 


36.64% 


0.92764 


0.88677 


18.53% 


1 


17 


0.05814 


0.09153 


57.75% 


0.10907 


0.08864 


43.24% 


1 


18 


0.04969 


0.08279 


30.12% 


0.06831 


0.0812 


17.81% 


1 


19 


0.10725 


0.08509 


52.91% 


0.08653 


0.07897 


19.23% 


1 


20 


0.0893 


0.08236 


62.2% 


0.06944 


0.08108 


21.85% 


1 


21 


0.02114 


0.01322 


47.88% 


0.01951 


0.01459 


30.43% 


1 


22 


0.01281 


0.01556 


50.13% 


0.01248 


0.01424 


25.29% 


1 


23 


0.01032 


0.01462 


31.32% 


0.01293 


0.01341 


18.77% 


1 


24 


0.0058 


0.01343 


38.54% 


0.00996 


0.0133 


20.48% 


3 


1 


4.01468 


3.04237 


34.49% 


2.93766 


3.06196 


16.15% 


3 


2 


2.43925 


2.88584 


33.77% 


2.67336 


2.9815 


12.16% 


3 


3 


1.66845 


1.96384 


36.29% 


1.94701 


2.05924 


24.65% 


3 


4 


1.61369 


1.89573 


45.73% 


1.79944 


2.00395 


17.04% 


3 


5 


0.39934 


0.17245 


46.26% 


0.22462 


0.17462 


19.1% 


3 


6 


0.30854 


0.15244 


34.22% 


0.26823 


0.17591 


19.79% 


3 


7 


0.14593 


0.17199 


35.86% 


0.17599 


0.1761 


19.23% 


3 


8 


0.3268 


0.15174 


46.17% 


0.22091 


0.17559 


29.68% 


3 


9 


0.03699 


0.02606 


45.78% 


0.03159 


0.02761 


21.52% 


3 


10 


0.0174 


0.02805 


46.25% 


0.0248 


0.02895 


22.52% 


3 


11 


0.03086 


0.02792 


47.53% 


0.0313 


0.02822 


21.98% 


3 


12 


0.01679 


0.02996 


35.62% 


0.02146 


0.02889 


15.67% 


3 


13 


0.92607 


1.08528 


113.38% 


0.88461 


1.01339 


44.1% 


3 


14 


0.66095 


0.86246 


41.28% 


0.89759 


0.93943 


25.98% 


3 


15 


0.32995 


0.57864 


39.28% 


0.71209 


0.63163 


17.53% 


3 


16 


0.76557 


0.59353 


66.16% 


0.80199 


0.63174 


27.11% 


3 


17 


0.10155 


0.05401 


50.01% 


0.05136 


0.05413 


33.84% 


3 


18 


0.08667 


0.05428 


85.39% 


0.06688 


0.05492 


35.01% 


3 


19 


0.05978 


0.0426 


57.18% 


0.05109 


0.04891 


26.68% 


3 


20 


0.03033 


0.05615 


121.29% 


0.06657 


0.0519 


34.31% 


3 


21 


0.00573 


0.01152 


102.81% 


0.00803 


0.01021 


44.1% 


3 


22 


0.00939 


0.00933 


100.7% 


0.01451 


0.01001 


43.49% 


3 


23 


0.01156 


0.00858 


102.21% 


0.00945 


0.00872 


37% 


3 


24 


0.00558 


0.0119 


83.98% 


0.00905 


0.00922 


34.99% 
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Tabic 9: Standard Monte Carlo and importance sampling estimates of sector contributions to VaR 
at 99.9% level and coefficients of variation of estimates. 



Sector 


Standard MC 


Importance sampling 




1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Cocf. of var. 


1 


13.09 


19.275 


23.61% 


21.336 


19.777 


7.92% 


2 


20.507 


20.512 


16.39% 


18.51 


19.683 


10.29% 


3 


17.853 


14.645 


23.53% 


13.302 


13.764 


14.01% 


4 


17.845 


14.784 


28.72% 


15.257 


15.474 


15.31% 



Table 10: Standard Monte Carlo and importance sampling estimates of portfolio-wide diversifica- 
tion index with respect to VaR at 99.9% level, in asset-level and sector-level context, and coefficients 
of variation of estimates. 



Context 


Standard MC 


Importance sampling 




1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Cocf. of var. 


Assets 


7.65% 


7.64% 


2.06% 


7.54% 


7.58% 


1.34% 


Sectors 


40.02% 


39.51% 


2.09% 


39.58% 


39.21% 


2.03% 
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Table 11: Standard Monte Carlo and importance sampling estimates of marginal diversification indices at 
asset-level with respect to portfolio-wide VaR at 99. 9% level and coefficients of variation of estimates. Only 
sectors 1 and 3. 



Sector 


Asset 


Standard MC 


Importance sampling 






1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Coef. of var. 


1 


1 


14.19% 


15.8% 


42.12% 


13.61% 


15.79% 


30.35% 


1 


2 


9.07% 


15.19% 


45.6% 


17.62% 


15.69% 


15.5% 


1 


3 


8.29% 


12.98% 


48.03% 


9.66% 


12.56% 


18.33% 


1 


4 


7.31% 


13.99% 


45.01% 


18.74% 


13.05% 


26.05% 


1 


5 


-0.85% 


5.66% 


75.78% 


4.97% 


4.27% 


31.38% 


1 


6 


4.15% 


4.34% 


84.72% 


6.48% 


4.67% 


25.68% 


1 


7 


-0.96% 


6.23% 


73.27% 


5.94% 


4.73% 


24.63% 


1 


8 


1.53% 


4.99% 


77.09% 


5.62% 


4.88% 


25.47% 


1 


9 


3.1% 


3.14% 


85.41% 


4.39% 


3.42% 


35.3% 


1 


10 


2.59% 


2.74% 


104.24% 


3.98% 


3% 


20.21% 


1 


11 


2.56% 


4.01% 


82.79% 


4.43% 


4.01% 


21.77% 


1 


12 


6.24% 


6.15% 


80.74% 


4.61% 


4.24% 


29.35% 


1 


13 


3.72% 


5.35% 


103.64% 


5.76% 


5.66% 


31.61% 


1 


14 


2.93% 


4.14% 


78.11% 


8.44% 


5.39% 


23.6% 


1 


15 


-0.38% 


3.15% 


95.93% 


4.61% 


5.53% 


30.56% 


1 


16 


4.98% 


4.37% 


71.17% 


6.87% 


4.96% 


21.72% 


1 


17 


0.12% 


1.86% 


145.91% 


1.67% 


1.81% 


67.75% 


1 


18 


0.28% 


2.07% 


114.54% 


0.93% 


1.53% 


28.81% 


1 


19 


-0.29% 


2.28% 


127.99% 


2.65% 


2.13% 


40.63% 


1 


20 


5.39% 


2.82% 


136.3% 


2.01% 


2.05% 


47.16% 


1 


21 


1.6% 


2.17% 


107.04% 


2.05% 


1.25% 


38.24% 


1 


22 


-0.21% 


0.84% 


243.25% 


1.18% 


1.37% 


50.75% 


1 


23 


1.48% 


2.12% 


135.81% 


1.45% 


1.79% 


26.88% 


1 


24 


-0.36% 


1.08% 


167.71% 


0.64% 


1.67% 


31.62% 


3 


1 


16.16% 


10.68% 


55.57% 


12.03% 


10.91% 


26.02% 


3 


2 


15.93% 


13.5% 


38.53% 


8.84% 


10.48% 


24.7% 


3 


3 


8.2% 


8.05% 


70.34% 


7.04% 


8.65% 


30.96% 


3 


4 


8.26% 


7.5% 


51.71% 


7.3% 


8.54% 


33.08% 


3 


5 


-0.9% 


2.03% 


129.42% 


5.16% 


2.42% 


44.94% 


3 


6 


4.44% 


3.83% 


88.67% 


5.06% 


2.19% 


44.09% 


3 


7 


3.35% 


2.94% 


108.77% 


2.17% 


3.04% 


39.1% 


3 


8 


6.85% 


2.28% 


118.93% 


4.82% 


2.71% 


52.58% 


3 


9 


2.65% 


1.86% 


141.32% 


2.08% 


1.63% 


49.99% 


3 


10 


3.19% 


3.01% 


102.98% 


1.15% 


1.8% 


52.75% 


3 


11 


0.33% 


2.21% 


148.43% 


2.77% 


2.15% 


56.96% 


3 


12 


1.27% 


2.2% 


135.55% 


1.2% 


2.47% 


36.1% 


3 


13 


5.88% 


5.16% 


92.98% 


3.78% 


4.37% 


91.65% 


3 


14 


7.2% 


4.31% 


82.07% 


3.31% 


3.59% 


31.49% 


3 


15 


3.55% 


2.94% 


105.14% 


2.9% 


3.18% 


33.27% 


3 


16 


6.08% 


4.54% 


81.95% 


4.11% 


3.42% 


47.75% 


3 


17 


-0.28% 


0.74% 


278.06% 


1.3% 


0.88% 


61.29% 


3 


18 


-0.27% 


0.79% 


241.35% 


1.12% 


0.9% 


81.37% 


3 


19 


-0.21% 


1.49% 


158.46% 


1.41% 


1.04% 


51.06% 


3 


20 


-0.38% 


0.71% 


254.37% 


1.07% 


1.13% 


108.81% 


3 


21 


2.31% 


1.5% 


157.1% 


0.27% 


0.9% 


89.92% 


3 


22 


5.26% 


0.31% 


435.57% 


1.17% 


0.79% 


83.06% 


3 


23 


-0.38% 


1.15% 


256.44% 


1.09% 


0.88% 


94.5% 


3 


24 


-0.38% 


0.61% 


248.39% 


0.74% 


1.17% 


78.71% 
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Tabic 12: Standard Monte Carlo and importance sampling estimates of diversification indices 
at sector-level with respect to portfolio-wide VaR at 99.9% level and coefficients of variation of 
estimates. 



Sector 


Standard MC 


Importance sampling 




1st run 


Mean of all runs 


Coef. of var. 


1st run 


Mean of all runs 


Cocf. of var. 


1 


30.46% 


44.51% 


25.27% 


51.01% 


45.94% 


9.39% 


2 


46.52% 


47.44% 


17.04% 


44.16% 


45.43% 


10.78% 


3 


41.37% 


32.77% 


25.98% 


30.24% 


30.76% 


15.13% 


4 


41.15% 


33.3% 


32.32% 


33.32% 


34.73% 


15.97% 
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